General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



X-763-70-353 

PREPRINT 


NASA m X- hS3Sf 


TIME AND FUEL OPTIMAL CONTROL FOR 
GRAVITY GRADIENT SPACECRAFT 


F. C. ZACH 


SEPTEMBER 1970 




GODDARD SPACE FLIGHT CENTER 

GREENBELT, MARYLAND 

IN 70 "4 20 I 8 


CN 

o 

■o 

2 

CL 

O 


(ACCESSION NUMBER) 


>SIONN 


(THRU) 


I 


(PAGES) _ 


i M_U 5 y ^L- 

u (NASA CR OR TMX OR AD NUMBER) 


(CODE) 

o I 


(CATEGORY) 


X-763-70-353 


PREPRINT 


TIME AND FUEL OPTIMAL CONTROL FOR 
GRAVITY GRADIENT SPACECRAFT 


F. C. Zach 


September 1970 


GODDARD SPACE FLIGHT CENTER 


Greenbelt, Maryland 


PRFced <nc page bun* k 

K N'OT F\Lmeo 

f mlvuuiI'.u rr.oL olai^n hoi r tuMcU. 


TIME AND FUEL OPTIMAL CONTROL FOR 


GRAVITY GRADIENT SPACECRAFT* 
F. C. Zach 


ABSTRACT 

Severe constraints on the amount of fuel available in spacecraft make it nec- 
essary to develop economical policies for attitude control. A promising approach 
is to combine active control by means of electric propulsion devices and passive 
gravity gradient control. The Maximum Principle is used to derive a time-fuel 
optimal control law for gravity gradient satellites. Switching lines which make 
it possible to determine the optimal control for each measured set of angles and 
rates are derived for decoupled gravity gradient satellites. It is shown how by 
this technique the special cases of pure time optimal and pure fuel optimal con- 
trol, which have been already solved in the literature, are covered. The possible 
extensions to coupled gravity gradient satellites are treated and it is shown that 
no unique switching lines can be derived. 


* The work for this paper was accomplished while the author held a National Research Council 
Postdoctoral Resident Research Associateship supported by the National Aeronautics and Space 
Administration, Goddard Space Flight Center. The author is very much obliged to Mr. W. Isley and 
Mr. D. Endres for their valuable suggestions. 
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NOMENCLATURE 

The index i stands for yaw, roll and pitch and is omitted when no distinction 
in the treatment of the different axes is necessary. 
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= gravity gradient coefficients, j = x, y, z 

= centers of the switching circular arcs, rr. = 1, 2, 3, 4; n = 1, 2, • ♦ • 
= ’ educed disturbance corque T di /I u ), i = x,y, z 
= Hamiltonian 

= moments of inertia, i = x , y , z 
= maximum absolute value of control u 
= performance index 
= auxiliary variables, m = 0, 1, • • • n 
= control torque, i = x,y, z 
= disturbance torque, i = x,y, z 
= time 

- time interval where control is on 
= time of control end 

= time interval where control is off 
= time of control start 

- reduced control torque (= T ci /I n ), i = x,y, z 
= state variable 

= cartesian coordinates of the centers oi the switching circular arcs, 
m - 1, 2, 3, 4; n = 1, 2, • • • 
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y m - y -coordinates of the switching circular arcs, m = 1, 2, 3, 4 
a = pitch angle (corresponds to z-axis) 



3 = maximum angle of circular arcs in thrusting periods 
r + = line for switching from u = 0 to u = + k 

P_ = line for switching from u = 0 to u = - k 

y = maximum angle of circular arcs in coasting periods 
A + = line for switching from u = k to u = 0 

A_ = line for switching from u = - k to u = 0 

K f = weighting factor for fuel in the j>erformance index P 
X t = weighting factor for time in the performance index P 
/c- , <j' m = polar coordinates of the centers of the switching arcs, m = 1, 2, 3, 4; 

n = 1, 2, • • • 

- roll angle (corresponds toy -axis) 

^ = yaw angle (corresponds to x-axis) 
co Q = orbital angular rate for circular equatorial orbit 




TIME AND FUEL OPTIMAL CONTROL FOR 


GRAVITY GRADIENT SPACECRAFT 


INTRODUCTION 

Attitude control of spacecraft requires economic use of fuei in order to 

maximize the number of control maneuvers for a specified amount of fuel. A 

promising approach is to combine active control by means of electric propulsion 

1,2 

devices and passive gravity gradient 'OG) control. Tt has been shown in Ref- 
erence 2 that the Maximum Principle can be used to derive the exact time optimal 
control ? i' / for gravity gradient satellites (GGS) if there is no coupling between 
the axes. However, if fuel constraints have to be imposed, consideration of fuel 
consumption in the optimum control law’ is necessary. An approximate solution 
has been derived in Reference 3 for GGSs with decoupled axes and in Reference 4 
for general GGSs. The approximations made in these references make it nec- 
essary to shift the switching lines in the phase plane according to the initial 
angles and rates. 

In the work performed here the exact solution for time/fuel optimal control 
of decoupled GGSs will be derived first. It will then be shown how the special 
cases of pure time optimal control as derived in References 2, 5, 6 and 7 and pure 
fuel optimal control (Reference 3) result from this approach. 

The possible extension of this time-fuel optimal control law to coupled GGSs 
will be given. 
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Control maneuvers treated here bring initial angles and rates which can be 
caused by disturbances to zero. Pointing in other directions can be treated as 
in Reference 2. 

TIME-FUEL OPTIMAL CONTROL FOR GGSs WITH UNCOUPLED AXES 
It has been shown in References 2 and 8 that the small angle equations of 
motion for a gravity gradient satellite in a nearly circular equatorial orbit with 
negligible inner damping can be written as 
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If it is assumed that 


^ i z ^ xx yy 


(4) 


the cross-coupling between che yaw and roll axes is eliminated. In this case the 
equations of motion can be written as: 


^ + a /, 'P = u ,/, + ^ + a <* * " u 4> + 
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It is assumed that the u . are constrained in magnitude, j.e. 


” k i £ u i £ i for i - a, 0, 0 
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If a disturbance free case(d - d f - d a - 0) and constant control (u , u f , 

u = constant) i3 considered, circles for the phase plots are obtained in the 

• • • _ 2 
(0, 0a 1 2 ), (£, 0a 7 1 2 ) and (a, aa 1 2 ) phase planes/ providing that a , a^ and 


a 2 t 0. The centers of these circles are at ik^a, with i = 0,0* a. For 




or a a =0 parabolas are obtained. For this case the solution for time-fuel 
optimal control has been derived in Reference 1. 

Equations (5) and (6) can be written as 


x + ax - u 


( 12 ) 


where x = 0,0, a, a = a 0 , a^ t a a , u = u^, u^, u a and“k<u<+k 


0 » 0 ! 
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5-7 

The Maximum lYinciple* will be applied to derive the time-fuel optimal 
control for Equation (12) in order to drive x and x to zero. The following pro- 
cedure has to be applied: 

(a) the performance index P has to be defined; for time-fuel optimal control 

P - I" (\ t +\ f |u|) dt ( 13 ) 

o 

where t is the start and t the end of control action. K x is the time weighting 
factor, \ f the fuel weighting factor. E.g. \ f = 0 leads to a pure time optimal 
solution. 

(b) the mathematical model of the system has to be written in the form of 
first order li.iear differential equations 

X, = f , (x. u) (i = 0, 1. • • • n) (14) 

n is the order of the system equations, x is the vector of the state variables 
(x Q • • • x n ) and u is the vector of the control variables (u Q • • *u n ). 
ic) the Hamiltonian 


H 



f i (x, u) 


i=0 

has to be maximized where p () • • • p n are auxiliary variables given by 
dp/dt = - 3H / d\ i and dx /dt = <9H/<9p i (i = 0, 1, ••• n) 


(15) 


(16) 
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It is shown in References 5—7 that for this case p 0 = - 1 can be chosen 
and *hat 


f 0 (*, u) = K % ♦ K f u 


( 17 ) 


In the case treated here 



— — 


n, — 


x o 


• X 

0 

% 

X 1 

\ 

x i 


_ X 2_ 




(18) 


which leads to 


K - p, x 2 + p 2 (u - axj ) - \ t - K f I u 


(19) 


H is maximized at every instant of time by 


u - ksgnp 2 for p 2 > ^ f 


( 20 ) 


and 


u - 0 for p 2 < 


( 21 ) 


With Equations (16) and (19) one obtains 


P 2 = P m Sin I 3 ' ^ 1 + ?) 


( 22 ) 


where p and cp are integration constants. 


Th e fina l arc: Figure 1 shows p 2 and u versus time t. From this figure it 
can be seen that u has to be on and off periodically. Furthermore, it can be said 
that u has to be on at the final time t e ; chis is true since the state would never 
reach the origin with control off because in this case the trajectories are circles 
with the origin as center. The on time before reaching the origin is between 0 
and t , t being the maximum on time. Therefore the part of the trajectory 
which leads the state to the origin has to be part of a semicircle with its center 
at 4 k/a or - k/aas shown in Figure 2. The longest arc leading to the origin is 
given by /3 = gj t c . 

Ten stingy before final arc : From Figure 1 it is obvious that the state has to 
have a coasting period before reaching this final arc unless the initial condition 
happened to be on this final arc. The length of the coasting period is limited by 
y - at jff . If the initial conditions are far enough from the origin that more thrust 
intervals than the final one leading to the origin are required to drive the state to 
the origin the coasting interval is given by 7 = ^t off . This means that by going 
back in time by t the time of switching from control to coasting has been 
found. In the phase plane this can be depicted by drawing circular arcs with the 
center ir .he origin and starting at each point of the arcs of Figure 2. The length 
of thes . arcs for coasting is given by y. The endpoints of these arcs again 
form circle arcs which are produced by rotating the circle arcs of Figure 2 by 
an angle counter clockwise about the origin. This is shown in Figure 3. It can 
be seen from this figure that each initial state lying within the hatched area can 
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be brought to the origin in a time-fuel optimal maneuver by coasting to \ or i 
respectively and then by control application. 

Thrusting before f inal coasting period: If the initial state does not lie in 
the hatched area of Figure 3 one sees from Figures 1 and 3 that another thrusting 
period is required wit.i control of the reverse sign when compared with the final 
control. This thrusting period again can last up to a time interval t c . Again 
shorter intervals than t for thrustings are applicable for cases where the initial 
condition is located such that it can be brought to by thrusting intervals less 
than t c . In the general case of thrust application of exactly t duration all pos- 
sible trajectories for this period again are circular arcs with an angle P and 
their centers at - k a or + k a, respectively. This means that the geometric locus 
of their starting points can be constructed by rotation of \_ or A + , respectively, 
in Figure 3 counterclockwise by an angle 3. The result is shown in Figure 4. 

Complete switching lines: The procedures outlined can be pursued in the 
same manner as above. The geometric relationship produced by rotating the 
circle arcs by or > , respectively, show' that all arcs which form the switching 
lines have their center on straight lines as shown in Figure 5. Also, the inter- 
section points of the arcs lie on a straight line. 

With y xa " 1 2 the centers C, of the circular arcs for n = 1, 2, • • • are 
given by 
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where 


In 


R([n(l + cos y) “ 1 ] 2 + [n sin y] 2 J 


(24) 




In 


a ret an In sin >/[n( 1 + cos > )- 1 ]} 


(25) 


R = k a 


(26) 


The equations of the circular arcs are then 


, = y,„ ' [ Rj -(*-* 1 .,)]' 


(27) 


for (n - 1 ) R( 1 + cosy) 1x1 nR( 1 + cos >) 


(28) 


To find the n corresponding to a given x the following equation can be applied 


n r Integer [x R(l + cosy)] + 1 (x>0) (29) 


The equations for the centers C 2n in the upper left hand quadrant are: 


X 2 n 

^2n COS ^2n ’ V 2 n 

■ Pin S1,1C P 2 n 

(30) 

■ 2 n 

r|[o( 1 + cos y) ~ cos y] 2 + [( 

(n " 1 ) sin y] 2 j 

(31) 

*2n 

= - arctan l(n - 1 ) sin')^[n( 1 

+ cos y) - cos y ] J 

(32) 

y c2 

= y 2 n + [r 2 - (x - X 2-1 ) 1 2 


(33) 
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for 


- nR( 1 + cos >)£x£-(n-l)R(l+ cos >) 


n for a given x can be found by 


Integer [ - x/R(l + cos y) + l] (x£0) 


The equations for the centers C 3n in the lower left hand ouaj^rant are: 


P 3n C °^3n * 


" sin>p 


^3n = *ln 


y 3 n + [r 2 - (x - X 3o ) 2 !' 


- nR( 1 + cos y) < x < - ( n - 1 ) R( 1 + cos y) 


n for a certain x is given by Equation (35). 


For the centers C 4n in the lower right hand quadrant the following equations 


hold 


- P c °scp 


4n ’ 


Ydn = Pin Sln <P, 
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( 41 ) 


4 n ^ 2n * 

y^ = y <n • '( x ' x 4„) 2 ]' 2 (42) 

for x given by inequality (28). 

n tor a certain x is given by Equation (29). 

Based upon Equations (23—42) the time/fuel optimal control law can be 
given by the flow chart of Figure 6. 

REDUCTION OF THE GENERAL TIME-FUEL OPTIMAL CONTROL TO PURE 
TIME AND PURE FUEL OPTIMAL CONTROL 

For the case where \ f = 0 the pure time optimal control law is derived. 

This means for the switching lines that y = 0 which leads to the well known 
semicircles as switching lines for the system wdth Equation (12). (Figure 7 and 
References 2, 5—7). 

On the other hand, for \ t = 0 the pure fuel optimal control law r is derived. 

In this case thrusting is only provided when the state is exactly on the xa~' 2 - 
axis, which means an impulse at this time. This corresponds wdth the results 
derived in Reference 3, where it is shown that most economic control for a sys- 
tem given by Equation 12 is achieved when the thruster is only in operation when 
the state falls on the xa' 1 2 - axis. 

OPTIMAL CONTROL OF GRAVITY GRADIENT SATELLITES WITH COUPLED AXES 
In the case w r here one principal moment of inertia is not equal to the sum of 
the two others the full GGS Equations (1)— (3) have to be considered. With 
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Equations (7)— (9) these can be written in the form: 


4> + % " ( l - *v/% 2 ) <t> = 

0 + 0 + (1 - a^4a 0 2 ) cc 0 <// = 


and 


a + a a - u + d 
a a a 


( 43 ) 


(44) 


(45) 


For the pitch equation the optimal control law derived above is still applicable; 
however, roll and yaw are now coupled. In order to calculate an optimal con- 
trol law the following transformations have to be made to make application of 
Equations (15) and (16) possible: 

a i = a j,/“'o 2 • a 2 = a V“o 2 (46) 

The time base shall be changed to r = cu Q t. Therefore 

• • • • 

d/dt = cu Q d/dr or 0 = 60 0 0 and 0 = 0 

which gives for Equations (43—44) 


and 


0 + a ! 0 + b t 0 ~ 


• • * — 

0 + a 2 0 + b 2 0 = 


(48) 


( 49 ) 
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with 


u 


u i/ lJ o 2 ( 1 = <*>) 


(50) 


and 


bj - - 1 , 


b 2 = 1 ' » 2 / 4 


(51) 


The following substitutions are made: 


0 = Xj , 


^ = x 2 , ^ - x 3 , 


0 - X 


(52) 


This yields for Equations (48—49): 


x i " X 2 


x 2 - - Xj - bj x 4 + u 0 


x 3 ' X 4 


X 4 = ' a 2 x 3 " b 2 X 2 + U <*> 


The performance index for time-fuel optimal control is 


(53) 


(54) 


P 


f [*-. * K f O^l + l u *l)] 


dt 


(55) 


which gives 


x 0 - K t + s d u J + i u ^i) 


(56) 
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The Hamiltonian is: 


H = " K t " K f (l u *l * l u j) + Pi x 2 * P 2 (" a i x i " b i x 4) 



+ P 2 + P 3 

X 4 + P* (" 

a 2 X 3' b J X 2 

) + P 4 u * (57) 

which is maximized by: 






sgnp 2 

for p 2 

2 A, 

(58) 


0 

for p 2 i 

< A, 

(59) 


k * s S n P 4 

for p 4 , 

> 

(60) 

% 

0 

for p 2 i 

< A, 

(61) 

Applying Equation (16) to Equation (37) yields: 



Pi ® 1 P 2 • P 2 

’ Pi + b 2 P4 ’ 

p 3 - a 

2 P 4 * P 4 

= _ P 3 4 b i P 2 (62) 

and 

•••• •• 
P 2 + A i P 2 

+ a 2 = 0 


(63) 


• ••• •• 

P 4 + A, p 4 

+ a 2 = 0 


(64) 


- 


V v 


with 


Aj - aj + a 2 - bj b 2 and A 2 = a t a 2 


(65) 


The roll and yaw equations can be written in a similar form if u. and u are 


constant: 


0 + A, 0 + A 2 - Oj 


( 66 ) 


0 + Aj 0 + A 2 - a 2 l V 


(67) 


The solutions for Equations (63) and (64) are: 


P 2 " P 21 sin (r + (p t ) + P 22 sin (A 2 j 2 r + qp 2 ) 


( 68 ) 


P4 = - ( b i/ b 2 ) P 2 l cos ( T+ '*>,) - (ai/a 2 )‘ 2 p 22 cos (A, 1 2 r + <P 2 ) (g9) 

. J 

For evaluation of the inequalities (58)— (61) only P 2] /P 22 (instead of P 21 and P 22 
separately) is interesting since \ f also is only interesting as its ratio to K for 
the performance index. To find trajectories similar to the case with decoupled 
axes a transformation is made for Equations (66) and (67) by 


Xj = $ + 0 , X 3 “ 'Jj' + 0 

■ j 

which results in 

*1 + A 2 X 1 = B 1 • *3 + A 2 X 3 = a 2 Q J. 



r 4 

■••4 

'.-4 


(70) 


(71) 
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Furthermore letting 


x i x 2 * X 3 X 4 

Equations (71) give circles in the (x Jf x 2 ) and (x,, x 4 ) phase planes for and u 
constant. 

To find the switching lines in these phase planes in analogy to the decoupled 
case, variation of <p J# <p 2 and P 21 /P, , in Equations (68) and (69) and fullfillment 
of the inequalities (58)— (61) as function of time in connection with the system 
time history in the (x { , x 2 ) and (x 3 , x 4 ) phase planes is required. 

However, no unique solution can be found. This i3 explained as follows: 

For the decoupled case the solution of the adjoint equation for roll and yaw 
is determined only by a phase angle cp (see Equation 22) which corresponds with 
the phase angles <p and <p in Equations (68)— (69). Due to coupling effects also 
P 21 /F 22 has to be given in order to evaluate inequalities (58)— (61). But P 21 /P 22 
not only depends on the coefficients bj and b 2 which determine the amount of 
crosscoupling but also on the initial conditions. The latter statement can easily 
be understood by the fact, that the switching sequence for one axis has to be much 
different whether the initial conditions on the other axis are large or small. This 
can be seen especially in these cases where the initial conditions are zero on 
one axis and ^ 0 on the other axis. In the decoupled case no control is required 
on the axis with initial conditions = 0. However, in the coupled case consider- 
able control action may be necessary also on the axis with zero initial conditions. 
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This fact has been the subject of experimental investigations on an analog 

computer and the optimum switching sequence for special cases has been found 

4 

by trial and error methods. 

CONCLUSIONS 

Time-fuel optimal control procedures have been derived for gravity gradient 
satellites where yaw and roll axes are decoupled. It has been shown that for the 
coupled case no switching lines can be found. (Switching lines are necessary for 
the practical application of an optimal control law because the control strategy 
lias to be based on the state of the system.) Therefore, it is recommended to 
try to decouple the yaw and roll axes of gravity gradient satellites by proper 
choice of moments of inertia if optimal control laws are to be applied. As one 
sees from the flow chart the control law is easy to apply. Choice of one parameter 
in the performance index covers the whole range between pure time and pure 
fuel optimal control. 
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1. Control u und auxiliary variable for time-fuel optimal control of the plant 
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Figure 2. Final Arc: This last part of the optin. . trajectory leads 
to the origin. The maximum duration of this arc is given 

by /2 - cut . 


COASTING 

REGION 



Figure 3. Coasting regions and switching lines. Between positive and negative 
thrust periods coasting periods are present whose length is given by y. 

+ • * • line for switching from u = 0 to + k 
’ _ • • • line for switching from u = 0 to ~k 


k + • • * line for switching from u = + k to 0 
' I ine for switching from u = ~ ’< to 0 
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Figure 6. Flowchart for implementation of the time fuel optimal control law 
INT (x) - integer of the number x. 
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